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Abstract: A knowledge of the tree-ring stable nitrogen isotope ratio (6!°N) can deepen our understanding 
of forest ecosystem dynamics by indicating the long-term availability, cycling and sources of nitrogen (N). 
However, the radial mobility of N blurs the interannual variations in the long-term N records. Previous 
studies of the chemical extraction of tree rings before analysis had produced inconsistent results and it is 
still unclear whether it is necessary to pre-treat wood samples from specific tree species to remove soluble 
N compounds before determining the 6'5N values. We compared the effects of pre-treatment with organic 
solvents and hot ultrapure water on the N concentration and 8!°N of tree rings from endemic Qinghai 
spruce (Picea crassifolia) growing in the interior of the central Qilian Mountains, China, during the last 60 a. 
We assessed the effects of different preparation protocols on the removal of the labile N compounds and 
investigated the need to pre-treat wood samples before determining the 5'4N values of tree rings. 
Increasing trends of the tree-ring N concentration were consistently observed in both the extracted and 
unextracted wood samples. The total N removed by extraction with organic solvents was about 17.60%, 
with a significantly higher amount in the sapwood section (P<0.01). The 8'5N values of tree rings decreased 
consistently from 1960 to 2019 in both the extracted and unextracted wood samples. Extraction with 
organic solvents increased the 15N values markedly by about 5.2%o and reduced the variations in the 6'!°N 
series. However, extraction with hot ultrapure water had little effect, with only a slight decrease in the 6!°N 
values of about 0.5%. Our results showed that the radial pattern in the inter-ring movement of N in 
Qinghai spruce was not minimized by extraction with either organic solvents or hot ultrapure water. It is 
unnecessary to conduct hot ultrapure water extraction for the wood samples from Qinghai spruce because 
of its negligible effect on the removal of the labile N. The 8'°N variation trend of tree rings in the 
unextracted wood samples was not influenced by the heartwood—sapwood transition zone. We suggest that 
the 8N values of the unextracted wood samples of the climate-sensitive Qinghai spruce could be used to 
explore the ecophysiological dynamics while focusing on the long-term variations. 
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1 Introduction 


As a promising environmental proxy, the tree-ring stable nitrogen isotope ratio (6!°N) has been 
used extensively to investigate long-term variations in the availability and cycling of nitrogen (N) 
in ecosystems (Robinson, 2001; Pardo and Nadelhoffer, 2010; Doucet et al., 2011; Gerhart and 
McLauchlan, 2014; Gessler et al., 2014; van der Sleen et al., 2017; Savard et al., 2020; Mason et 
al., 2022). However, the reliability of using the 5!°N series of tree rings has been impeded because 
of the potential influence of the mobility of the unstable N components within plant stems and the 
low amounts of N stored in the wood tissue (Cutter and Guyette, 1993; Meerts, 2002). The radial 
movement of the unstable N fractions has restricted the accuracy of the interannual resolution of 
the tree-ring N series. The N translocation mechanism is a response to the internal 
physiological processes that sustain tree growth. N is carried from the older tree rings to the 
cambium tissue during the lignification process with the senescence of the sapwood (Meerts, 
2002). With the death of xylem parenchyma cells, N in cytoplasm is retrieved for possible reuse 
elsewhere in the trees as the sapwood is transformed into the heartwood (Merrill and Cowling, 
1966). The movement of unstable N masks the actual N isotope signal in tree rings formed in a 
given year and the precise historical information about soil N dynamics (Mizota et al., 2011). 
Therefore, wood samples are often extracted before analysis to try to overcome this problem. 

Wood extraction experiments, which were aimed at removing the labile N components before 
tracing the actual N signature of the annual tree rings, have yielded contradictory results among 
different tree species. Sheppard and Thompson (2000) reported that extraction with organic 
solvents reduced the amount and variability of the total N concentration of tree rings from 
ponderosa pine (Pinus ponderosa) and concluded that this method successfully removed mobile N 
compounds. However, these results were not replicated when the same extraction protocol was 
applied to different tree species. For example, extraction with organic solvents removed the 
unstable N in the sapwood, but did not affect the overall trend of the N concentration in tree rings 
from red pine (Pinus densiflora) (Kwak et al., 2009). The N concentration and 5!°N values in tree 
rings from Norway spruce (Picea abies) were not significantly affected when they were extracted 
using a chemical protocol (Tomlinson et al., 2014). 

A recent study concluded that solvent-based extraction did not effectively eliminate the 
non-structural N in tree rings of western redcedar (Thuja plicata) because the !°N-labeled signal 
was detected in tree rings formed up to about 10 a previously (Bunn et al., 2017). The influence of 
the wood extraction experiment on the same species has been shown to be inconsistent; for 
example, solvent-based pre-treatment of wood increased the interannual resolution of the N 
concentration and ô! N values, but did not remove all mobile N in beech (Fagus sylvatica), 
however, in another investigation, pre-treatment of wood decreased the variation in the 5!°N 
signal (Elhani et al., 2003, 2005). Many researchers have modified the original 
Sheppard—Thompson protocol, such as simplifying the original method, changing the organic 
solvents and using different pre-treatment methods (Hart and Classen, 2003; Bukata and Kyser, 
2005; Savard et al., 2009; Doucet et al., 2011; Drake et al., 2011; Larry et al., 2011). In summary, 
there is still considerable uncertainty in the assessment of wood pre-treatment methods for certain 
tree species used to record potential climate and environmental signals. 

In this study, we documented the N concentration and 5!°N values of tree rings in endemic 
Qinghai spruce (Picea crassifolia), one of the dominant coniferous tree species growing on shady 
and semi-shady slopes in the Qilian Mountains, China (Liang et al., 2006). We used three 
different pre-treatment methods before measuring the N concentration and 6'°N values of tree 
rings: (1) wood extracted with organic solvents; (2) wood extracted with hot ultrapure water; and 
(3) unextracted wood. The aim of this study was to investigate the ability of different extraction 
methods to remove the soluble N-containing compounds from wood samples. We explored 
whether it is necessary to pre-treat wood samples to give reliable 5!°N series for tree rings of 
Qinghai spruce for use in investigating the historical record of N availability and cycling in 
temperate forests in remote mountainous areas where N deposition is low. 
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2 Materials and methods 


2.1 Study site and field sampling 


The study site is located at the Sidalong Forest Farm (99°54'10"E, 38°26'58"N) in the upper 
reaches of the Heihe River catchment in the interior of the central Qilian Mountains, China. The 
climate in this region is influenced by both the East Asian Monsoon and the Westerlies (Li and 
Liu, 2000). Climate data from two nearby meteorological stations (Qilian and Minle) showed that 
the annual mean temperature in the time period of 1958—2019 ranged from 1.0°C to 4.1°C and the 
total annual precipitation varied from 288.2 to 526.6 mm. 

We collected two cores per living tree from a forest stand of pure Qinghai spruce at an 
elevation of 2590 m a.s.l., close to the lower timberline of Qinghai spruce. We totally collected 54 
tree ring cores from opposite sides of the stem of 27 large mature Qinghai spruce trees in June 
2020. The samples were taken at about 1.3 m above the ground (breast height) using a 10-mm 
diameter increment borer. 


2.2 Wood pre-treatment 


2.2.1 Sample preparation 


We collected all the tree cores to the laboratory and processed them according to the standard 
dendrochronology techniques (Stokes and Smiley, 1996). After air-drying and polishing, we used 
the LINTAB 6.0 system (RINNTECH, Heidelberg, Germany) for ring width measurements and 
applied the COFECHA program (Holmes, 1983) to guarantee the quality of cross-dating. One 
core from each of five trees was chosen for manual dissection into annual rings using a clean 
stainless-steel blade with the aid of a binocular microscope (SDPTOP, Suzhou, China). We pooled 
the annual rings from five cores of the same year into one sample for isotope analysis. The mixed 
samples were ground in a ball-mill to a powder that could pass through a 200-mesh screen. All the 
powdered samples were oven-dried to constant weight at 55°C. 


2.2.2 Nitrogen (N) extraction protocol 


Each annual powder sample was separated into three subsamples (weight of 55—80 mg for each): 
(1) wood extracted with organic solvents; (2) wood extracted with hot ultrapure water; and (3) 
unextracted wood (Fig. 1). 

Unstable compounds (e.g., resins, salts and most free forms of alkaloids and tannins) were 
removed from the subsamples in the first group via solvent-based extraction following the 
protocol of Sheppard and Thompson (2000). The subsamples were loaded into filter bags and 
extracted for 10 h in a 1:2 mixture of toluene and ethanol and then for 10 h in ethanol in a Soxhlet 
apparatus. The filter bags were then rinsed three times for 2 h each time with ultrapure water (18 
MQ.cm) until the resultant solution became transparent. 

The second group of subsamples was extracted with hot ultrapure water to remove soluble 
compounds, including tannins, inorganic salts and low molecular weight polysaccharides (Rowell, 
1984; Ajuong and Breese, 1997) following a modified version of the water-bath method of Bukata 
and Kyser (2005). We used hot water rather than cold water because some compounds (e.g., 
tannins) are more soluble at higher temperatures. A 1 mL volume of ultrapure water (18 MQ.cm) 
was added to each centrifuge tube containing the subsamples and heated at 90°C for 4 h in a dry 
block bath (THB-2, AS ONE Corp., Osaka, Japan). The tubes were centrifuged at 10,000 r/min 
for 5 min and the supernatant solution was decanted after standing for 1 h. 

All the extracted and unextracted subsamples were then oven-dried at 55°C for 48 h to remove 
any remaining water and reagents. 


2.3 Determination of the stable N isotopes 


Wood subsamples of about 7 mg were enclosed in tin capsules and combusted in an Element 
Analyzer (EA Isolink, Thermo Fisher, Bremen, Germany) coupled via a ConFlo VI interface to a 
Gas Isotope Ratio Mass Spectrometry (Delta V Advantage, Thermo Fisher, Bremen, Germany) at 
the Laboratory of Stable Isotope and Global Change in Shaanxi Normal University, Xi'an, China. 
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Fig. 1 Flow chart of the wood pre-treatment methods for the determination of the stable nitrogen isotope ratio 
(SN) in tree rings from the endemic Qinghai spruce. Three pre-treatment methods were used to prepare the wood 
samples for the measurements of the nitrogen (N) concentration and 5!°N in tree rings: (1) wood extracted with 
organic solvents (solvent-extracted wood); (2) wood extracted with hot ultrapure water (water-extracted wood); 
and (3) unextracted wood (bulk wood). 


We used the certified international standard material USGS-55 (8!°N= —0.3%o (£0.4%o)), which 
has a high N concentration (0.25%), for instrument calibration. To eliminate any nonlinear effects 
associated with the small peak heights, we developed a wood internal standard collected from 
Qilian juniper (Sabina przewalskii) trees, which could produce peak heights similar to those from 
the international standard. We calibrated the international standard using material USGS-55 to 
assess the accuracy of our method (standard deviation lower than 0.3%0; n=128). The stable N 
isotope ratio was expressed as per mil (%o) value using delta (ò) notation as relative deviations 
from the international standard (atmospheric N2), and was normalized using one-point linear 
interpolation. 

Because the N concentration in tree rings was very low relative to the carbon (C) concentration, 
we used an ascarite trap (12 mL volume) to scrub excessive CO2 from the sample gas before gas 
chromatography. The analysis time was set at 600 s to eliminate potential memory effects 
originating from the large amount of CO2 that was produced from preceding combustion. 


2.4 Statistical analyses 


The heartwood—sapwood boundary is not clearly detectable by differences in the color of the 
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wood from Qinghai spruce. The most abrupt shift in the N concentration of tree rings during the 
time period of 1960-2019 was determined using the changepoint package in the R statistical 
program (version 4.1.3). We detected an abrupt increase in the N concentration in the late 1990s, 
close to rings from 1995 to 1998, suggesting that the wood of this species might contain a 
transition zone of the heartwood to the sapwood in the physiological sense of the term (Fig. 2) 
(Merrill and Cowling, 1966). 

To minimize the differences and evaluate the similarities between the three curves for the 
different pre-treatment methods, we normalized the absolute values of the N concentration and 
SDN in each series to obtain the standard Z-scores. Differences in the N concentration or 5'°N 
values among the three different pre-treatments were determined by one-way ANOVA followed 
by a least significant difference (LSD) test of means with a 95% confidence level. We also 
checked some differences in the series of the N concentration and 5'°N between the sapwood and 
heartwood and/or among different extractions by exploiting some statistical indices. Differences 
in the offsets of the N concentration or 5!°N values between the extracted wood and unextracted 
wood were compared by an independent sample t-test. Differences in the N concentration or 5!°N 
values between the sapwood and heartwood among the three different pre-treatments were also 
tested via an independent sample t-test. Statistical analyses were completed with SPSS Statistics 
25 (IBM, USA) and performed in R 4.1.3 (R Core Team, 2022). 
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Fig. 2 N concentration (a) and 6'°N values (b) in tree rings from Qinghai spruce as a function of time 
determined by the three different pre-treatment methods: wood extracted with organic solvents (solvent-extracted 
wood), wood extracted with hot ultrapure water (water-extracted wood) and unextracted wood (bulk wood). The 
vertical gray-shaded area indicates the heartwood—sapwood transition zone around the time period of 1995-1998. 


3 Results 


3.1 Tree-ring N concentration 


Increasing trends of the tree-ring N concentration were observed in Qinghai spruce for the 
extracted wood (solvent-extracted wood and water-extracted wood) and unextracted wood (bulk 
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wood) from 1960 to 2019. There was a typical distribution pattern in which the N concentration 
of tree rings increased sharply in the heartwood—sapwood transition zone, whereas there was a 
constant flat trend in the heartwood, followed by a steady increase in the sapwood, with the 
highest concentration in the most recent rings (Fig. 2a). The year-to-year variation in the N 
concentration of the sapwood was greater than the variation in the N concentration of the 
heartwood, in which the interannual correlations were stronger, as observed in some of the 
statistical indices (Table S1). The variations in the N concentration of the solvent-extracted wood 
and water-extracted wood were consistent with those in the bulk wood, as indicated by the 
significant linear correlations between the standardized Z-scores of the extracted wood and 
unextracted wood (R?=0.955 between the solvent-extracted wood and bulk wood, and R2=0.969 
between the water-extracted wood and bulk wood; P<0.01) (Fig. 3). 
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Fig. 3 Comparisons of the time series of the tree-ring N concentration and 5!°N values for the extracted wood 
(solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood) during the time period of 
1960-2019, as well as relationships of the Z-scores of the N concentration and 6'°N values between the extracted 
wood (solvent-extracted and water-extracted wood) and unextracted wood (bulk wood). (a and b), the N 
concentration; (c and d), the 5!5N values. The curves in the left panel were standardized by the Z-scores. The solid 
lines in the right panel represent the linear regression of each dataset. The 1:1 line is shown by the dashed line. 
The coefficient of determination (R2) and P-value are also given. 


Pre-treatment of wood with organic solvents removed about 17.60% of the total N, with a 
statistically significant difference between the heartwood and sapwood sections (P<0.01). Tree 
rings of the solvent-extracted wood had remarkable lower N concentration (mean of 0.052%, 
P<0.05) than tree rings of the water-extracted wood and bulk wood (means of 0.066% and 
0.063%, respectively), as shown by the results of one-way ANOVA followed by an LSD test 
(Table 1; Fig. 4). The differences in the N concentration between the solvent-extracted wood and 
bulk wood were significantly larger in the sapwood section than in the heartwood section 
(P<0.01), with mean values of —0.016% and —0.008%, respectively, as indicated by the 
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independent sample t-test (Table 1; Fig. 5). The coefficient of variance reflected the lower 
interannual variability in the N concentration of the solvent-extracted wood compared to the bulk 
wood (Table S1). 

Extraction by hot ultrapure water had little effect on the removal of the soluble N compounds 
from tree rings. There was no significant difference in the N concentration in tree rings between 
the water-extracted wood and bulk wood (P>0.05; Fig. 4). Differences in the N concentration 
between the whole wood extracted by ultrapure water and bulk wood increased marginally with a 
mean value of 0.002% (Table 1). No statistically significant difference was observed between the 
heartwood and sapwood in the water-extracted wood and bulk wood from the results of the 
independent sample t-test (P=0.127; Table 1; Fig. 5). 


Table 1 Summary statistical analyses among the pre-treatment methods of wood samples 


N concentration (%) 


Bulk wood Solvent-extracted wood = Water-extracted wood Df. I Df. H 
Whole wood 0.063+0.017 0.052+0.012 0.066+0.017 —0.011+0.006 0.002+0.003 
Sapwood 0.081+0.017 0.065+0.012 0.084+0.014 —0.016+0.005 0.003+0.004 
Heartwood 0.053+0.008 0.045+0.004 0.055+0.007 —0.008+0.005 0.002+0.002 
P 0.000°* 0.000°* 0.000™* 0.000** 0.127 
SIN (%o) 
Bulk wood Solvent-extracted wood = Water-extracted wood Df. I Df. H 
Whole wood 8.4+4.1 13.6+3.7 7.8+3.9 5.2+2.1 —0.5+0.8 
Sapwood 4.141.8 9.52.4 3.71.5 5.5+0.8 —0.3+0.6 
Heartwood 11.142.8 15.9+1.7 10.5+2.6 4.942.6 —0.6+0.9 
P 0.000°* 0.000°* 0.000™* 0.322 0.228 


Note: Mean+standard deviation values of the tree-ring N concentration and ôN values for the extracted wood 
(solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood). Df. I and Df. II refer to the difference 
between the solvent-extracted wood and bulk wood and between the water-extracted wood and bulk wood, respectively. The P 
value represents result from the independent sample t-test, which is dimensionless. **, P<0.01 level. 
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Fig. 4 Box—whisker plots of the tree-ring N concentration (a) and 6!°N values (b), as well as relationships 
between the N concentration and 6!5N values (c) of the extracted wood (solvent-extracted wood and 
water-extracted wood) and unextracted wood (bulk wood). The boundaries of the boxes indicate the 25" and 75" 
percentiles. The solid lines are median values, the black dots are mean values and the black squares are outliers. 
Bars represent the change range of the N concentration or 5'°N values. Different lowercase letters indicate 
significant differences (P<0.05) between the pre-treatment methods by one-way ANOVA. The solid curves 
represent the correlations simulated by exponential fitting. 
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Fig. 5 Box—whisker plots of the changes in the tree-ring N concentration (a) and 5'°N values (b) for the 
heartwood and sapwood sections in the extracted wood (solvent-extracted wood and water-extracted wood) 
compared with the unextracted wood (bulk wood), as well as relationships between differences in the N 
concentration and 8!°N values between the extracted wood and unextracted wood (c). The boundaries of the boxes 
indicate the 25'* and 75" percentiles. The solid lines are median values, the black dots are mean values and the 
black squares are outliers. The bars represent the change range of differences in the N concentration or 6°N 
values. Offsets were calculated between the extracted wood and unextracted wood. Different lowercase letters 
indicate significant differences (P<0.05) between the sapwood and heartwood by the independent sample t-test. 


3.2 Tree-ring ðN 


There was a consistent decrease during the research period in the 5!°N of tree rings from Qinghai 
spruce for all three pre-treatment methods (Fig. 2b). Temporal variations in the tree-ring 5'°N of 
the solvent-extracted wood and water-extracted wood corresponded with those in the bulk wood, 
indicated by significant linear relationships between the normalized Z-score curves of the 
extracted and unextracted wood (R?=0.734 for the solvent-extracted wood and R?=0.965 for the 
water-extracted wood; P<0.01; Fig. 3). The results were similar to those for the N concentration: 
the variations in the tree-ring 5'°N in the sapwood were larger than those in the heartwood and the 
interannual correlations were stronger in the heartwood (Table S1). However, we also observed 
some differences in the tree-ring 6!°N signal among the three pre-treatment methods. 

The removal of the mobile N by extraction with organic solvents weakened the variations in the 
tree-ring 5!°N signal, but markedly increased the mean (5.2%o) of the absolute values (Table 1). 
The mean tree-ring 5!°N value of the solvent-extracted wood (mean of 13.6%) was mostly higher 
than those of the water-extracted wood and bulk wood (7.8%o and 8.4%o, respectively), which 
were significantly different (P<0.05), as indicated by the ANOVA LSD test (Table 1; Fig. 4). 
However, there was a gentle trend in the tree-ring 5'°N of the solvent-extracted wood compared 
with the bulk wood, as indicated by the slopes of linear regression lines (Fig. $1; Table S2) and 
the values of the coefficient of variance (Table $1). The independent sample t-test showed no 
significant difference (P=0.322) in the offsets of the tree-ring 6!°N values between the heartwood 
and sapwood for tree rings extracted with organic solvents (Table 1; Fig. 5). These results 
mirrored those of the solvent-extracted wood, in which the tree-ring 5'°N values in the sapwood 
did not increase significantly relative to those in the heartwood. 

Pre-treatment with hot ultrapure water had little effect, with the mean tree-ring 5!°N value 
decreasing slightly by 0.5%o. The 5'°N values of the water-extracted wood decreased marginally 
and insignificantly with a mean of 7.8% (P>0.05), as shown by the ANOVA LSD test (Table 1; 
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Fig. 4). Although differences in the tree-ring 5'°N values between the water-extracted wood and 
bulk wood were higher in the heartwood than in the sapwood (mean values of —0.6%o and —0.3%o, 
respectively), no significant difference (P=0.228) was observed by the independent sample f-test 
(Table 1; Fig. 5). 


4 Discussion 


4.1 Higher efficiency of chemical extraction in removing mobile N 


Given the efficiency of removing the mobile N compounds, extraction with organic solvents may 
be more popular than extraction with hot ultrapure water. Solvent-based extraction reduced the 
total N concentration in the whole wood by about 17.60% and removed significantly more N in 
the sapwood section of tree rings, whereas the offsets in the N concentration between the 
water-extracted wood and bulk wood were almost zero (Table 1; Fig. 5a). Many researchers have 
reported that the tree-ring N concentration is lower after solvent-based extraction following the 
Sheppard—Thompson protocol to remove the unstable N (Sheppard and Thompson, 2000; Elhani 
et al., 2003, 2005; Kwak et al., 2009; Hietz et al., 2010; Tomlinson et al., 2014), although the 
proportion of the N removed differs among tree species (Fig. 6). For instance, the total N 
concentration of tree rings in the sapwood decreased significantly by about 45.00% for ponderosa 
pine (P. ponderosa) (Sheppard and Thompson, 2000) and decreased from 0.22% to 0.15% after 
extraction for beech (F. sylvatica) (Elhani et al., 2003). For red pine (P. densiflora), extraction 
with organic solvents removed 19.40%-—31.60% of the total N in recent tree rings (Kwak et al., 
2009), in agreement with our results that solvent-based extraction significantly reduced the total 
N in the sapwood section of tree rings from Qinghai spruce (P. crassifolia) by about 19.30%. 
Differences in the N concentration between the water-extracted wood and bulk wood were close 
to zero in both the heartwood and sapwood (Table 1), contrasting with the finding that the weight 
percentage of the N decreased significantly by 10.50% in the sapwood of tree rings extracted with 


m= Sapwood mE Heartwood mm Whole wood 


Qinghai spruce Qinghai spruce 
White oak White oak 
Japanese black pine Japanese black pine 


Big-leaf mahogany Big-leaf mahogany 


Spanish cedar Spanish cedar 


Red pine Red pine 


Beech Beech 


Douglas-fir Douglas-fir 


Ponderosa pine Ponderosa pine 


50 —40 -30 -20 -10 0 10 20 


Offset in the 3'5N (%o) i poe aa 


Fig. 6 Offset in the tree-ring 5'°N values (a) and percentage of offset in the N concentration (b) for different tree 
species after the extraction of bulk wood. Offset in the 8!°N was calculated as the difference between the 
extracted wood and unextracted wood, so offset in the N concentration did. The percentage of offset in the N 
concentration was calculated as the quotients of the offset in the N concentration divided by the actual value of 
the unextracted wood. Tree rings in white oak were extracted with deionized water. Tree rings of other species 
were extracted with organic solvents. Data were obtained from previous studies as follows (see Table S3 for 
details): [1], Sheppard and Thompson (2000); [2], Elhani et al. (2003); [3], Elhani et al. (2005); [4], Kwak et al. 
(2009); [5] Hietz et al. (2010); [6], Larry et al. (2011); [7], Bukata and Kyser (2005); [8], Hart and Classen 
(2003). S-e, solvent-extracted wood; W-e, water-extracted wood. 
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deionized water in non-resinous species such as white oak (Quercus alba) (Bukata and Kyser, 
2005). The complex influences may be caused by the diversity of the wood extractives in tree 
rings from different tree species. Resins form the principal part of the wood extractives of 
temperate conifers (Mutton, 1962; Fengel and Wegener, 1989; Sheppard and Thompson, 2000) 
and are not extracted with hot water as a result of their hydrophobic nature. 


4.2 Distinct differences in the tree-ring 5'°N values between the two extraction methods 


Pre-treatment of the bulk wood samples with organic solvents reduced the variations in the 
tree-ring 5!°N signal (Fig. S1; Table S2) and remarkably increased the absolute values (Table 1), 
in line with the results of the non-labeled beech trees (F sylvatica) extracted by organic solvents 
(Elhani et al., 2005). The higher 5!°N values for tree rings in the extracted wood have also been 
reported for Spanish cedar (Cedrela odorata), big-leaf mahogany (Swietenia macrophylla), 
Japanese black pine (Pinus thunbergii) and Norway spruce (P. abies) (Hietz et al., 2010; Larry et 
al., 2011; Tomlinson et al., 2014). Extraction with organic solvents increases the absolute values 
of the 6!°N in the wood samples, which could be explained by assuming that the decomposition of 
the N compounds into mobile forms depletes the '°N of the soluble N compounds as a result of 
fractionation (Yoneyama et al., 1998). Solvent-based extraction fractionates the °N so that the 
lighter '4N is removed, which could explain the elevated 5'°N values of tree rings after extraction 
(Tomlinson et al., 2014). 

The 5'°N values of tree rings extracted with ultrapure water as a function of time was almost 
coincident with those for the bulk wood (Fig. S1; Table S2), although the absolute values 
decreased slightly (Table 1; Fig. 4). Our results are consistent with previous reports that there are 
no significant changes in the tree-ring 5!°N values for red oak (Quercus rubra) and white oak 
(Quercus alba) after pre-treatment in water-bath (Bukata and Kyser, 2005). Extraction with hot 
ultrapure water had little effect on the 5'°N signal of tree rings from Qinghai spruce and did not 
change the 5!°N values beyond the limit of analytical precision (£0.3%o). 


4.3 Proposal of using the bulk wood for the determination of the 6!5N in tree rings of 
Qinghai spruce 

Extraction with either organic solvents or hot ultrapure water did not eliminate the 
physiologically driven radial pattern in the N concentration associated with the N translocation, 
suggesting that not all of the mobile N fractions were removed from tree rings. Although 
solvent-based extraction significantly reduced the total N concentration in the wood, the 
distribution pattern of the N concentration was not minimized. The tree-ring N concentration and 
5°N values for the solvent-extracted wood were probably underestimated because extraction with 
organic solvents not only removed the soluble N-containing compounds, but also other 
compounds. An unexpected finding was that the translocation of N between inter-rings had a 
visible impact on the 5!°N series of tree rings in the solvent-extracted wood. The major decrease 
in the 5'°N values over the whole series coincided with a sharp increase in the N concentration in 
the heartwood—sapwood transition zone (Fig. 2a). The effect of water extraction on either the N 
concentration or 6!°N values was negligible. Hence, neither extraction with organic solvents nor 
hot ultrapure water eliminated the N translocation pattern, which might mask the original annual 
SN signal of tree rings. 

The 5°N values of the bulk wood in Qinghai spruce should be a more reliable method of 
determining long-term variations in forest ecosystems. The most obvious variation in the 6!°N 
series of tree rings in the bulk wood did not occur in the heartwood—sapwood transition zone (Fig. 
2), suggesting that the physiological pattern in the N concentration does not influence the 6'°N 
values. This could provide indirect evidence that the radial mobility of N among tree rings does 
not generate significant fractionation of N isotopes in Qinghai spruce (Doucet et al., 2011, 2012). 
Given the confusing effect of different pre-treatments in the removal of the mobile N compounds, 
we recommend that the bulk wood is used for the determination of the 6!°N in tree rings from 
endemic Qinghai spruce to provide long-term environmental signals. 
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4.4 Consistent downward trend of the tree-ring 64°N in different protocols 


Variations in the 5!°N of tree rings from all three pre-treatment methods showed a consistent 
decrease for Qinghai spruce trees growing in the interior of the central Qilian Mountains during 
the research period (Figs. 2b and 3). The level of the N deposition at our study site was very low. 
The average atmospheric dry and wet N deposition values over the region were 9.6 and 9.5 kg 
N/(hm?-a) during the time period of 2011-2015, respectively, both lower than the average fluxes 
across China (20.6 and 19.3 kg N/(hm?.a), respectively) (Xu et al., 2015; Jia et al., 2016, 2019). 
As a result of the limitations in the temporal resolution of N deposition data, we cannot detect 
reliable direct or indirect evidence of anthropogenic atmospheric N deposition in the annual 
variations of the 5'°N in tree rings, in contrast with the 6!°N values of tree rings from riparian 
trees in the lower reaches of the Heihe River, China (Wang et al., 2020). The more pronounced 
decreasing trend of the 6!°N in tree rings might be caused by fractionation in the soil—forest N 
cycle after the forest was hold allochthonous N inputs with the lower 6!°N, which requires further 
investigation. We ruled out the possibility that the downward trend in the tree-ring 5!°N in recent 
decades might be related to different soil textures or soil depths (Högberg, 1997) because there 
were no significant local differences in soils affecting the sampled trees from this pure, mature 
Qinghai spruce stand. 

The decreased trend of the 5'°N values in tree rings showed that the availability of N in the 
ecosystem had decreased and that the N cycle in these soils has become increasing closing over 
the last 60 a (Austin and Vitousek, 1998; Emmett et al., 1998; Bukata and Kyser, 2005; 
McLauchlan et al., 2007; Craine et al., 2015; McLauchlan et al., 2017). Tree growth, as indicated 
by the basal area increment of Qinghai spruce, has accelerated since the 1950s in the interior of 
the central Qilian Mountains as a result of a warmer and wetter climate (Figs. S2 and S3). A 
longer growing season resulting from a warmer climate and increased atmospheric CO2 
concentrations could stimulate a greater N demand for tree growth, which could reduce the 
availability of N (Ollinger et al., 2002; Norby et al., 2016). Mineralization, which is the dominant 
source of the bioavailable N for trees growing in temperate forests in remote mountain regions, 
generally increases under wet conditions (Guerrieri et al., 2010; Greaver et al., 2016). Ecosystems 
might therefore meet higher plant demands for N by decreasing the N losses and increasing the N 
fixation and mineralization, leading to a gradually closing N cycle (Peñuelas and Estiarte, 1997). 

Unfortunately, we do not have evidence of long-term trends to examine whether our series of 
the N concentration and 6'°N in tree rings could be related to the tree ages. The increases in the 
carbon isotope ratio (6!°C) and/or reductions in the intrinsic water use efficiency within the first 
few decades of the rapid growth of trees most likely reflect the age effects, which are greatly 
influenced by physiological changes linked to the height growth during tree development 
(Peñuelas et al., 2008; Härdtle et al., 2013; Guerrieri et al., 2020). We speculated that variations 
in the N concentration and 5!°N of tree rings during their juvenile period might be more dramatic, 
similar to the large interannual variability in the sapwood of mature trees. Living xylem 
parenchyma cells sustain rapid growth in young trees before the lignification process is complete, 
affecting physiological processes within the tree stem (Merrill and Cowling, 1966). 

There were nonlinear relationships between the N concentration and 6!°N values of tree rings in 
Qinghai spruce for all three pre-treatment methods (Fig. 4c). The negative correlation was 
stronger in the sapwood section than in the heartwood section, consistent with the discovery that 
this negative relationship decreased when excluding the sapwood samples (Trumper et al., 2020). 
These differences in the correlation between the heartwood and sapwood may help to explain the 
differences in this relationship under lower or higher N concentration. Nevertheless, we have 
insufficient evidence to demonstrate whether this relationship is meaningful. If a more 
appropriate extraction technique is developed to remove the N translocation mechanism, then the 
N concentration of tree rings could be used as another environmental indicator. The potential 
meaning of the correlation between the N concentration and 6'°N values in tree rings would then 
be of more interest to researchers. 
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5 Conclusions 


Extraction based on organic solvents may be more efficient in removing the mobile N 
components from the wood of endemic Qinghai spruce species than hot ultrapure water. Our two 
extraction experiments produced a distinct difference in the 5!°N of tree rings in both the 
tendency and absolute values. Not all the unstable N compounds were removed by any 
pre-treatment of the wood samples, which could be explained that the results of wood 
pre-treatment with organic solvents differ strongly among tree species. 

Given that the impact of wood pre-treatments on the tree-ring N isotopes is complex and varies 
between species, we recommend that the bulk wood samples of tree rings from Qinghai spruce 
are used to register long-term variations in the environment. To increase the temporal resolution 
and dependability of the 5'°N series in tree rings, we suggest: (1) the modification of current 
and/or development of new pre-treatment methods to improve the efficiency of extraction; (2) 
corroboration of whether the removal of the soluble N compounds is necessary in view of 
species-specific differences; and (3) the verification of evidence to directly demonstrate whether 
differences in physiological processes affect the environmental signals recorded in tree rings. 
Although we are still unable to pinpoint the exact mechanism underlying the decreased 5'°N in 
tree rings, a triple stable isotope approach combining 65°C, 6'8O0 and 65N_ with 
dendrochronological analyses has potential in monitoring long-term physiological responses of 
Qinghai spruce trees to global climate change. 
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Appendix 


Table S1 Summary statistical analyses, including mean sensitivity, 1‘ autocorrelation and coefficient of 
variance of the N concentration and stable nitrogen isotope ratio (8!5N) in tree rings of the extracted wood 
(solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood) 


N concentration 


Mean sensitivity 


1“ autocorrelation 


Coefficient of variance 


Wood section 


Bulk Solvent- Water- Bulk Solvent- Water- Bulk Solvent- Water- 
extracted extracted extracted extracted extracted extracted 
wood wood wood 
wood wood wood wood wood wood 
Whole wood 0.083 0.073 0.084 0.702 0.668 0.728 0.272 0.236 0.254 
Sapwood 0.101 0.091 0.111 0.376 0.325 0.262 0.208 0.188 0.169 
Heartwood 0.059 0.052 0.051 0.807 0.732 0.829 0.143 0.092 0.121 
"N 
Mean sensitivity 1“ autocorrelation Coefficient of variance 
Wood section 
Bulk Solvent- Water- Bulk Solvent- Water- Bulk Solvent- Water- 
extracted extracted extracted extracted extracted extracted 
wood wood wood 
wood wood wood wood wood wood 
Whole wood 0.139 0.115 0.135 0.876 0.765 0.889 0.494 0.271 0.492 
Sapwood 0.262 0.168 0.307 0.472 0.385 0.259 0.450 0.254 0.394 
Heartwood 0.101 0.088 0.086 0.828 0.452 0.874 0.256 0.108 0.244 
0.15 (a) — Solvent-extracted wood 0.15 5 (b) — Solvent-extracted wood 
— Water-extracted wood — Water-extracted wood 
So — Bulk wood — Bulk wood 
5 0.12 — -= -—- Regression line O12 ssaa Regression line (heartwood) 
g saute anaes omnes Regression line (sapwood) 
g 
= 0.09 0.09 
3 
f= 
8 
z, 0.06 0.06 
0.03 0.03 
1960 1980 2000 2020 1960 1980 2000 2020 
G 
w — Solvent-extracted wood ¥\: ; 
— Solvent-extracted wood —— Water-extracted wood 
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— Bulk wood ee Regression line (heartwood) 
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1960 1980 2000 2020 1960 1980 2000 2020 
Year Year 


Fig. S1 Trends simulated by the linear regression of the N concentration (a) and &!°N (c) in tree rings of the 
extracted wood (solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood), and 
trends of the N concentration (b) and 8'°N (d) of the heartwood and sapwood sections in tree rings of the extracted 
wood (solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood) from 1960 to 2019. 
The vertical grey-shaded area indicates the heartwood—sapwood transition zone around the period of 1995-1998. 
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Table S2 Statistical analyses of the variations in the N concentration and 6!°N in tree rings of the extracted 
wood (solvent-extracted wood and water-extracted wood) and unextracted wood (bulk wood) 


N concentration (%) 


Wood Heartwood Sapwood Whole wood 
Slope R Slope R? Slope R? 
Solvent-extracted wood 0.00030 0.532 0.00096 0.230 0.00056 0.629 
Water-extracted wood 0.00058 0.808 0.00077 0.110 0.00082 0.730 
Bulk wood 0.00068 0.831 0.00113 0.165 0.00082 0.689 
S!5N (%o) 
Wood Heartwood Sapwood Whole wood 
Slope R? Slope R Slope Re 
Solvent-extracted wood —0.08 0.216 —0.15 0.129 —0.17 0.659 
Water-extracted wood —0.23 0.853 —0.07 0.092 —0.21 0.888 
Bulk wood —0.25 0.802 —0.12 0.167 0.22 0.869 


Note: Variations in the heartwood, sapwood and whole wood are represented by the slopes of the regression lines (see Fig. S1). The 
coefficient of determination (R°) is provided. 


Table S3 Retrospect of extraction protocols used in the tree-ring 5!°N analysis from previous studies 


Method Extraction protocol Tree species N concentration 65N Determination Reference 
: The variation in N 
Ponderosz arda 
ce A Puls Decrease (about concentration of tree Sheppard and 
(Pinus ponderosa) D NA : ; Thompson 
$ na 45.000%) rings was reduced 
and Douglas-fir ; (2000) 
substantially. 
D ase (about The interannual resolution 
ecrease (abou = 
of the N and 5N 
Beech (Fagus 36.000%, from Decrease (about was improved, but not Elhani et al. 
sylvatica) nee 0.4%0) all mobile N Ja (2003) 
E removed. 
Increase (—6.0%o— Variation in the 6!°N signal 
—5.0%0 to was decreased, but eee 
S1:4hina1:1 Ee A NA —4.5%o -3.5%0 differences between trees a etal, 
apie Fy f REPES, in the control before and after treatment 
ee trees) were enhanced. 
Sheppard/ S2:4hin than Decrease Higher extractable N ii the 
£ ` ý Red pine (Pinus sapwood was removed, Kwak et al. 
Thompson : i (19.400%-— NA 
S3: 1 or 4 h in densiflora) B 600%) while not affecting the (2009) 
distilled i á overall trend of N. 
water Spanish cedar 
(Cedrela odorata) Decrease Īncrease (about The proportion of N A 
and big-leaf aie ‘ z extractable in the Hietz et al. 
(not significantly 1.0%o higher in i 
mahogany lowe the treated wood) sapwood was not higher (2010) 
(Swietenia than that in the heartwood. 
macrophylla) 
Norway spruce No'signiticant effect Tomlinson et 
(Picea abies) Decrease Increase (slight) a te tree-ring N and al. (2014) 
15 
Western redcedar Hievated. X The mobile N in tree rings 
ee signal found baat Bunn et al. 
(Thuja plicata) NA was eliminated 
and Douglas-fir 10 a before the ineffectivel: (2017) 
Ens 15N-labeled y 
Increase of 3.0%o 
Simplified Keep overnight ina Japanese black aes a Significant effect on the iaevieval 
Sheppard/ 1:1 mixture of pine (Pinus NA deea of tree-ring 5!°N under (20 it) : 
Thompson toluene/ethanol thunbergii) en oa high-N conditions. 
—7.0%bo in the sites 
with the high N 
S1:in 1:1 mixture of Beech (F. Decreas The difference is always 
Altered benzene/methanol; He difo iae (8.000% an d Decrease (0.4%o within the analytical error, Doücetetai 
Sheppard/ S2: 12 or 48 h in ai ae 7 pein and 0.3%o, with no significant change (2011) ` 
Thompson acetone; S3: 1 hin (Picea ssi s) respectiva) respectively) in the N or 65N values or 


deionized (DI) water 


temporal trends. 


To be continued 


202206.00162v1 


chinaXiv 


ChinaXiv@ (ERAT! 


WANG Ziyi et al.: Tree-ring 5'°N of Qinghai spruce in the central Qilian Mountains... 


Continued 

Methods Extraction protocol Tree species N concentration 8°N Determination Reference 

Fresh DI water Red oak (Quercus Decrease Increase ed 
Water-bath (shaking 4 times each rubra) and white  (10.500% in the (0.6%o in the niece change A 

for 3 d) oak (Quercus alba) sapwood) sapwood) ` y 

S1: 16-18 h ina 2:1 Decrease acne eae an 

: 4 jist 

Holocellulose mixture of toluene/ Ponderosa pine Decrease (2.0%0-3.0%o in control trees, while Hart and 
extraction ethanol; S2: 16-18 h the control trees; relatively laroeand Classen 

in ethanol; $3:6hin (ponderosa) (about 1.000%) 4 994, 34.0% in "C AUVElY arg (2003) 

DI water the labeled trees) Y2tlable effect on the 

highly !5N. 


Note: NA, not analyzed. 
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Fig. S2 Growth pattern of Qinghai Spruce during the period of 1840-2019. Variations in the tree growth were 
estimated as the basal area increment (BAI) during the period of 1840-2019. The orange line is smoothed using a 
20 a low-pass Fast Fourier Transform (FFT) filter to emphasize the low-frequency variations. 
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Fig. S3 Temporal trends of annual mean temperature (a) and total precipitation (b) from 1958 to 2019 based on 
averaging data from two nearby meteorological stations (Qilian and Minle) of the study area. The dashed lines 
represent the trends simulated by the linear regressions. Annual variation is represented by the slope of the trend 
line. The coefficient of determination (R°?) and significance P-value are also provided. 
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